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, ^ , Abstract. We present a calculation of thermal photon production i.e. photons from 

secondary interactions among particles produced in heavy ion collisions at collider 
energies. This is done within the framework of hydrodynamics. We take into account 
the lack of chemical equilibrium in QGP. It turns out that main effects from chemical 
non-equilibrium composition of QGP, reduction of particle number and increase in 
temperature, nearly cancel in photon spectrum. 
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2 ■ 1. Introduction 

^ I Thermal emission of real photons is very sensitive to the initial temperature of emitting 
matter and their mean free path in the produced matter is so long that they escape 
^ , from the fireball without further interactions. This makes thermal photons an excellent 
^ I probe of the early stage of the collision. Unfortunately, even in best situations, thermal 
^ ' photons are only a small contribution to the total photon spectrum. At low transverse 
momentum kx, the main contribution is from the electromagnetic decays of hadrons, 
^ I mostly from tt", and at high kr from prompt photons i.e. photons from primary 

\ interactions of the partons of the colliding nuclei. At collider energies, however, we 
have a possible window at intermediate /cr — 3 — 7 GeV to detect the thermal part 
of the photon spectrum. Also a new contribution to the photon spectrum has been 
proposed recently As high energy jets propagate through QGP they will induce 

emission of real photons. Photons from this process will mix with the thermal spectrum 
increasing the overall signal from QGP. 

One uncertainty in thermal emission is the chemical composition of the QGP. 
We present here a calculation of the thermal part of the photon spectrum within the 
framework of hydrodynamics Supplemented with photon emission rates in thermal 
matter, hydrodynamics provides a closed framework to study the thermal emission in 
a heavy-ion collision. Within this framework it is relatively easy to take the chemical 
composition of the QGP into account. It turns out that the effects from chemical 
non-equilibrium, the decrease of quark number relative to gluons and the increase in 
temperature, nearly cancel in the thermal photon contribution to the spectrum. 
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2. Theoretical framework 

Our framework to study the thermal photon emission is the following: The initial energy 
and particle densities are obtained from the pQCD + saturation model |1] which also 
gives the production time of the initial parton matter. Assuming thermal equilibrium 
immediately after production, the whole expansion of the matter can be treated by using 
hydrodynamics. This model, with the extra assumptions of longitudinal boost invariance 
and that the matter behaves as an ideal fluid, is in good agreement with the measured 
hadron spectra at RHIC Since we study the effects of particle composition on the 
photon emission, we relax the assumption of chemical equilibrium. We also neglect the 
small net baryon number. 

Chemical composition of the matter can be controlled by introducing chemical 
potentials which we approximate by using multiplicative fugacities for both quarks and 
gluons. The time evolution of fugacities is determined by rate equations and following 
Biro et al. P, processes qq gg and gg ^ ggg are included. Rate equations can 
be solved together with the hydrodynamic equations [El El EI, once the equation of 
state is specifled. This model gives a complete space-time evolution of energy density 
and transverse velocity as well as temperature, fugacities and other thermodynamic 
quantities. Once the space-time evolution is known from hydrodynamics and the photon 
emission rate as a function of temperature and fugacities is specifled, we can integrate 
the rate over the space-time history of the collision to obtain the spectrum of thermal 
photons. Photon emission rate in the hot hadron gas is calculated and parametrized 
in refs. [HI E3 E]- In QGP in chemical equilibrium the emission rate, complete to 
leading order in aem and ag, is calculated in [T2]. Following the method of ref. [TH], this 
calculation has been extended to the non-equilibrium case in ref |14j . 

3. Equation of state 

The equation of state (EoS) is specifled by combining the hadron gas equation of state 
and the QGP bag model EoS by Maxwell construction. The hadron gas is an ideal gas 
of all hadronic states up to a mass 1.3 GeV and QGP is treated as an ideal gas of gluons 
and three flavors of massless quarks. In the approximation of multiplicative fugacities 
Xi for quarks and gluons the particle densities in the QGP are given by 



where hi is the particle density in chemical equilibrium. Within the same approximation, 
pressure is given by 



where B is bag constant. The connection between the critical temperature and the 
bag constant B is given by Gibbs phase co-existence condition 




(1) 



(2) 




(3) 
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where PQGP,th is the thermal part of the pressure in QGP Q. 

The use of fugacities for describing the gluon and quark densities out of chemical 
equilibrium leads to an ambiguity on how to choose the bag constant and the critical 
temperature (cf. eq. (3)). In the full kinetic and chemical equilibrium with fugacities 
equal to one the situation is clear since the temperature is the only independent variable 
(neglecting net baryon number) in both phases. With QGP out of chemical equilibrium, 
we do not know quantities in the hadron gas which would be uniquely defined from the 
quark and gluon fugacities at phase transition. Here we fix the ambiguity in matching 
by simply taking Tc to be a constant, independent of fugacities. In other words, we 
assume that the QGP always hadronizes at the same temperature into the same state of 
hadron gas. In our calculations = 167 MeV both in and out-of chemical equilibrium. 
Because the photon emission from QGP is dominated by the early, high density stage 
of the collision our results are not sensitive to the details of the phase transition. 

4. Results 

The time evolution of temperature and fugacities in the center of the fireball are shown in 
the Fig. 1. As can be seen from the top figure the matter is close to chemical equilibrium 
at the end of the QGP phase. This justifies further our choice of fixed value for Tc. In 
the lower figure we can see the time evolution of temperature compared with the full 
chemical equilibrium case. Same initial energy density is used both in equilibrium and 
in out-of-equilibrium case. 

These two figures show the effect of lack of chemical equilibrium. When the energy 
density is kept constant, reduction of the particle number below its equilibrium value 
leads to an increase of temperature. Under these conditions there are less particles but 
they have higher average energy. These deviations from the equilibrium values have 
opposite effects on the photon emission rates and nearly cancel each other. This can 
be seen from Fig. 2, where we plot thermal photon spectrum for both cases. Only at 
high transverse momenta, ^ 10 GeV, the increase in temperature is becoming more 
important and the equilibrium spectrum starts to fall faster than the non-equilibrium 
one. However the experimentally accessible region is around A;^ ~ 5 GeV, where the 
spectrum is practically independent of the chemical composition of QGP. 

5. Discussion 

We have calculated the spectrum of thermal photons, emitted during the expansion stage 
of heavy ion collision, in the case when the QGP is not in chemical equilibrium. The 
space-time evolution of the heavy ion collision is modeled by hydrodynamics together 
with rate equations for the time evolution of chemical composition of QGP. Thermal 
photon spectrum is then obtained by integrating the temperature and in the QGP 
the fugacity dependent photon emission rates over this space-time history. We have 
demonstrated that the lack of chemical equilibrium does not change much the photon 
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Figure 1. Time evolution of fugacities and 
temperature in the middle of the firebah in the 
LHC lead-lead collision. 
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Figure 2. Thermal photon spectrum for the 
LHC lead-lead collision both in and out-of 
chemical equilibrium. 



spectrum as compared to the chemical equilibrium case. The reason for this is that 
while there are less particles in the non-equilibrium QGP, the temperature becomes 
higher. These have opposite effects on photon emission and nearly cancel in the region 
where the possibilities to detect experimentally the thermal part of the spectrum are the 
best. Thus in heavy ion collisions thermal photons are not measuring initial temperature 
alone, but the combination of temperature and particle densities. 
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